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ABSTRACT 

Modern (sub-)millimeter/radio interferometers such as ALMA, JVLA and the PdBI successor NOEMA will 
enable us to measure the dust and molecular gas emission from galaxies that have luminosities lower than 
the Milky Way, out to high redshifts and with unprecedented spatial resolution and sensitivity. This will pro- 
vide new constraints on the star formation properties and gas reservoir in galaxies throughout cosmic times 
through dedicated deep field campaigns targeting the CO/[CII] lines and dust continuum emission in the (sub- 
millimeter regime. In this paper, we present empirical predictions for such line and continuum deep fields. 
We base these predictions on the deepest available optical/near-infrared ACS and NICMOS data on the Hubble 
Ultra Deep Field (over an area of about 12 arcmin 2 ). Using a physically-motivated spectral energy distribution 
model, we fit the observed optical/near-infrared emission of 13,099 galaxies with redshifts up to z = 5, and 
obtain median likelihood estimates of their stellar mass, star formation rate, dust attenuation and dust luminos- 
ity. We combine the attenuated stellar spectra with a library of infrared emission models spanning a wide range 
of dust temperatures to derive statistical constraints on the dust emission in the infrared and (sub-)millimeter 
which are consistent with the observed optical/near-infrared emission in terms of energy balance. This allows 
us to estimate, for each galaxy, the (sub-)millimeter continuum flux densities in several ALMA, PdBI/NOEMA 
and JVLA bands. As a consistency check, we verify that the 850/im number counts and extragalactic back- 
ground light derived using our predictions are consistent with previous observations. Using empirical relations 
between the observed CO/[CII] line luminosities and the infrared luminosity of star-forming galaxies, we infer 
the luminosity of the CO(1-0) and [CII] lines from the estimated infrared luminosity of each galaxy in our 
sample. We then predict the luminosities of higher CO transition lines CO(2-l) to CO(7-6) based on two 
extreme gas excitation scenarios: quiescent (Milky Way) and starburst (M82). We use our predictions to dis- 
cuss possible deep field strategies with ALMA. The predictions presented in this study will serve as a direct 
benchmark for future deep field campaigns in the (sub-)millimeter regime. 

Keywords: dust, extinction - galaxies: evolution - galaxies: ISM - galaxies: statistics - submillimeter. 



1. INTRODUCTION 

The last decade has seen impressive advances in our under- 
standing of galaxy formation and evolution through galaxy 
surveys done (preferentially) at optical and infrared wave- 
lengths. In particular, the history of star formatio n (the 
'Lillv-Madau' plot: e.g. Lilly et alj (l99"6l: lMadau et alJ l!996h 
Hopkins & Beacom 2006), and the build up of stellar mass 
as a function of galaxy type and mass, have been well quan- 
tified to within 1 Gyr of the Big Bang. It has been shown 
that the comoving cosmic star formation rate density likely 
increases gradually from z ~ 6 — 10, it peaks at z ~ 2, 
and drops by more than an order of magnitude from z ~ 1 
to z ~ dHopkins & Beacomll2006l: iBouwens et alJl2010h . 
The build-up of stell ar mass follows this evolution, a s does 
the temporal integral dBell et al.112001 lllbert et aT1l2010l) . The 
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redshift range z ~ 1 — 3, during which roughly half of 
the stars in the Universe were formed, is referred to as the 
'epoch of galaxy assembly'. In summary, the star forma- 
tion properties as well as the stellar masses of galaxies have 
been well delineated through these optical and near-infrared 
deep field studies; almost all of our current knowledge is 
based on optical and near-infrared deep field s of the stars, 
star formation, and ioni zed gas (but see, e.g. Smolcic et al. 
2009, Kari m et alJl201 lL for additional constraints based on 
radio continuum studies). E.g., Lyman Break Galaxy samples 
have reve aled a major popula tion of star-forming galaxies at 
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Likewise, magnitude -selected 
12001 iLillv et al.ll2007h provide 



a census of the star-forming population based on UV/optical 
flux rather than color. 

A key measurement that is currently (mostly) unavailable 
is that of the presence of molecular gas, i.e. the dense ISM 
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phase ('fuel') required for star formation to take place, which 
lies at the heart of the evolution of the comoving cosmic 
star formation rate density. In recent years, there have been 
significant efforts devoted to obtaining molecular gas mea- 
surements of individual galaxies, by performing follow-up 
studies of galaxies that have been pre-selected from opti- 
cal/NIR deep surveys. To date, color-selection techniques 
(e.g., 'BzK', 'BMBX') have revealed significant samples of 
gas-rich, star forming galaxies at z ~ 1.5 to 2.5 (with molec- 
ular gas masses Mh 2 — 10 10 — 10 11 M Q , stellar masses 
M* ~ 10 10 - 10 11 M m and star for m ation rates SFR ~ 
100 M ? yr-^lDaddi et al.ll2008ll2010bSlGenzel et alJ lIoTol: 
iTacconi et al.ll2010tlGeach et al.ll201 ll) . While very important 
in their own right, these studies (that focus on the detection of 
carbon monoxide, the main tracer for molecular gas at low and 
high redshift) remain fundamentally limited to galaxy popu- 
lations that were pre-selected in the optical/near-infrared, i.e. 
potentially missing gas-dominated and/or obscured systems. 

From a_ theore tic al/mod eling perspective, 

lObreschkow et alj d2009alfbl 1201 lh provide simulations 
of the cosmic evolution of the molecular (and atomic) 
hydrogen in galaxies as a function of re dshift, by building on 
the Millennium dark matter simulations (Springel et al.ll2005h 
in which they place 'idealized model galaxies' at the centers 
of the dark matter halos which then ev olve according to 
simple r ules (' se mi-analytica l mode ling', [Obreschko w et al.l 
I2009allbl l20lll) . iPower et al.1 d20lob and lLagos et all (i20TTl) 
also present models of the cosmic evolution of the atomic 
and molecular gas content in galaxies by applying different 
semi-analytical galaxy formation models to the Millennium 
simulation. 

In this paper, we present empirical predictions of molecular 
line and continuum deep fields that are only now becoming 
possible thanks to the advent of observational facilities that 
dwarf previous capabilities, in particular the broad bandwidth 
and sensitivity afforded by the Atacama Large Millimeter Ar- 
ray (ALMA), the Jansky Very Large Array (JVLA, formerly 
known as EVLA) and the IRAM Plateau de Bure Interferome- 
ter (PdBI) successor, the Northern Extended Millimeter Array 
(NOEMA). Our predictions are based on the deepest available 
optical and near-infrared data available for the Hubble Ultra 
Deep Field (UDF), but, barring cosmic variance, should give 
a statistical representation of an arbitrarily chosen region on 
the sky. Basically, we use a sophisticated spectral energy dis- 
tribution (SEP) mod el combined with a Bayesian approach 
(Ida Cunha et al.ll2()08l) to interpret the observed optical/near- 
infrared emission of the UDF galaxies in terms of their stellar 
content, star formation activity and dust attenuation, and ob- 
tain statistical constraints on their total dust luminosity which 
are consistent with the observed stellar emission in terms of 
energy balance (i.e. all the stellar radiation absorbed by dust 
in the rest-frame ultraviolet to near-infrared is re-emitted in 
the mid-infrared to millimeter range). The dust luminosity 
of each model is then re-distributed at infrared to millimeter 
wavelengths by combining the (dust-attenuated) stellar SED 
with a library of infrared dust emission models spanning a 
wide range of dust SED shapes (including different dust tem- 
peratures). This allows us to derive median-likelihood esti- 
mates of the (sub-)millimeter continuum flux densities of our 
galaxies in several ALMA, PdBI/NOEMA and JVLA bands. 
Based on these continuum predictions, we calculate predicted 
line strengths in the various rotational transitions of carbon 
monoxide (CO) and ionized carbon ([CII]), two main tracers 
of the star-forming interstellar medium in galaxies. We note 



that, with this technique, we potentially miss, by definition, 
extremely dust obscured sources that are not included in our 
optical/near-infrared catalog. However, this should not have a 
great impact in our results since the main goal of this paper is 
to characterize the general galaxy population rather than the 
extreme, dust-enshrouded starbursts. 

In Section |2j we describe the main optical/near-infrared 
photometric catalogue of the Hubble UDF in which we 
base our predictions, as well as additional data from opti- 
cal, infrared and sub-millimeter surveys of the UDF/Extended 
Chandra Deep Field South area that we use to test our predic- 
tions. In Section [5J we describe our spectral energy model 
and fitting method. In Section[4] we analyze in detail our pre- 
dicted (sub-)millimeter properties for 13,099 galaxies in the 
UDF. We discuss the SED fitting outputs and derived p hys- 
ical prop erties of the galaxies in our sample in Sections |4. 1 1 
and 14.21 respectively. In Section l4~3l we discuss the reliabil- 
ity of our infrared luminosity and (sub-)millimeter continuum 
flux estimates fr om the obser ved optical/near-infrared SEDs, 
and in Sections 14.41 and 14.51 we perform consistency checks 
on the predicted continuum flux densities of our galaxies at 
850/im by comparing them with observed nu mber counts and 
the extragalactic background light. In Section l4~6] we present 
the distribution of our predicted continuum flux densities for 
the whole sample in sever al (s ub-)millimeter bands from 38 
to 870 GHz, and in SectiongZZlwe infer the CO and [CII] line 
luminosities of the galaxies in our sample from the estimates 
of their infrared luminosities, based on simple empirical re- 
lations. Based on these (sub-)millimeter line and continuum 
predictions, we discuss future deep fields with ALMA in Sec- 
tion [5] We summarize our main conclusions in Section [6] 
Appendix lAl contains a comparison of our results with what 
we would obtain assuming fixed SED shapes for the galaxies 
in our sample. 

Throughout this paper, we use a concordance ACDM cos- 
mology with Hq = 70 km s _1 Mpc -1 , S1a = 0.7 and 
fi m = 0.3. 

2. THE DATA 

As an example, we here base our predictions on the deepest 
available optical and infrared data on the Hubble Ultra Deep 
Field (UDF) that will be accessible with ALMA. We stress 
that this coherent database is the only reason for our choice 
and that, barring the issue of cosmic variance, we could have 
used any other field for our predictions as well. This means 
that our statistical predictions should also hold for a northern 
field that will be accessible from other telescopes (e.g. IRAM 
PdBI/NOEMA, JVLA). 

2.1. Input catalog: optical/near-infrared HST data 

We start by using the photometric catalog of the Hubble 
UDF (centere d at R.A.= 03 h 32 m 39 s .O, Dec=-27°47'29".l) 
described in ICoe et all (120061) . This contains aperture- 
matched, PSF-corrected ACS BVi'z' and NICMOS3 JH 
photometry, as well as Bayesian photometric redshifts for 
all the detected sources, accurate to within 0.04(1 + z spcc ) 
(ICoe et al.ll2006l) . The full catalog contains 18, 700 sources, 
of which a large fraction (8,042) are detected at the IOct level 
in at least one band. The IOct limiting AB magnitudes in the 
B, V, i', z', J and H are 28.71, 2 9 13, 29.01 28.43 , 28.30 
and 28.22, respectively. Following [Coe et ail d2006l) . to ex- 
clude contamination from stars, we exclude sources with i- 
band stellarity stel > 0.8 (about 6 per cent of the sample), 
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Figure 1. Distribution of the photometric redshift (a) and V-band magnitude 
(b) of the optically-selected Hubble UDF galaxies from the lCoe et aTfoOQg) 
catalog. 



leaving us with 17,532 galaxies, with photometric redshifts 
distributed as plotted in Fig. [TJ a )- F° r reference, in Fig. [TJb), 
we plot the distribution of the observed ACS V-band magni- 
tude for 16,830 sources detected in that band. We note that the 
fact that the redshift distribution of our sources peaks at z ~ 2 
and the sudden drop in sources fainter than 30 AB magnitudes 
in the V-band are due to incompleteness. This implies that 
our predictions may be missing high-redshift, dust-obscured 
galaxies that are too faint in the optical to be in our sample and 
may have large (sub-)millimeter fluxes. This is the case of 
LESSJ0333243. 6-274644, the only sub-mm source detected 
in the UDF as part of the LESS survey (LABOCA observa- 
tions of the Exte nded Chandra Deep Field South at 870 /im; 
IWeifi et aill2009l) . which has no optical counterpart in our op- 
tical catalog. Our predictions are therefore lower limits for 
the possible number of detections at high redshift (z > 2). 

2.2. Supporting data 

We complement the photometry in the UDF catalog with 
additional photometry out to the far-infrared. We use 54 
galaxies detected in the Herschel/SPIRE bands available 
in the publicly r eleased HerMES survey (P.I. S. Oliver; 
lOliveretal .1120101) images in GOODS-South, for which we 
applied the same prior source extraction technique as in 
lElbaz et al.l d2011l) for the GOODS-Herschel SPIRE data in 
GOODS-North. HerschelfPACS images of the GOODS- 
South are available as part of the GOODS-Herschel program 



(P.I. D. Elbaz). For each of these 250/im-selected galaxies, 
well-sampled spectral energy distributions from the ultravio- 
let to the far-infrared are available, including photometry in 
the following bands: U, B, V, i, z, J, K, SpitzerllRAC 3.6, 
4.5, 5.8 and 8.0 /im, Spizer/IRS at 16 /im, SpitzerfMTPS at 
24 pm, HerschelfPACS at 70, 100 and 160 /im, and Her- 
schel/SPIRE at 25 0, 350 and 500 pm (lElbaz et al.l I20TH 
Magdis et al. 201 1). The redshifts of the galaxies in this sub- 
sample go from z = 0.140 to z = 2.578, wit h a m edian of 
value of 1.0. We use this sub-sample in Section l4~3l to test the 
reliability of our predictions of the total infrared luminosity 
and (sub-)millimeter continuum fluxes from observed optical 
data as described in Section [3] 

To test our predictions for a wid er fie ld and address the is- 
sue of cosmic variance in Section 14.41 we use a wider-area 
catalog of the Chandra Deep Field South field which also in- 
cludes the UDF but is abo ut 10 time s larger in area, the FIRE- 
WORKS catalog fWuyts et alJl2008h . The FIREWORKS cat- 
alog is a A' s -band selected galaxy catalog which contains 
multi-wavelength photometry of 6,308 galaxies from the U 
band to the Spitzer-24pm band, with a 5cr magnitude limit of 
24.3 AB mag in the K s band (i. e. s hallower than the UDF 
HST catalog described in Section l2~TT >. over 138 arcmin 2 . 



3. SED MODELLING 

In a next step we use the models described in 
Ida Cunha et al.l (120081) to fit the observed rest-frame optical 
to near-infrared spectral energy distributions of the gala xies 
from the photometric catalog described in Section 12.11 and 
obtain statistical estimates of the infrared luminosities, (sub- 
millimeter continuum flux densities and CO line luminosities 
for each galaxy in the sample. 

3.1. Ultraviolet to near-infrared emission 

We use the spectral synthesis model of Bruzual & Charloj 
(120031) to compute the integrated light emitted by stars in 
galaxies for a wide range of metallicities (distributed between 
0.02 and 2 times solar), ages (distributed between 0.1 Gyr 
and the age of the Universe at each redshift), and star forma- 
tion histories (parameterized as exponentially declining with 
a wide range of timescales, and superimposed random bursts 
of star formation). To account for the attenuation of starlight 
by dust, we describe the interstell ar medium of galaxie s us- 
ing the two-component model of ICharlo t & Fall (2000|): the 
ambient (diffuse) interstellar medium and the star-forming re- 
gions (birth clouds). This prescription accounts for the fact 
that stars are born in dense molecular clouds, which dissi- 
pate typically on a time-scale of 10 7 yr. As a result, the non- 
ionizing continuum emission from young OB stars and line 
emission from their surrounding HII regions is absorbed by 
dust in these birth clouds and then in the ambient ISM, while 
the light emitted by stars older than 10 7 yr propagates only 
through the diffuse ISM. The main free parameters of this 
model are the effective V-band optical depth seen by young 
stars in birth clouds, fy, and the fraction of fy contributed 
by dust in the diffuse ISM, p. Using this model, we compute 
the attenuated stellar emission of galaxies and the total lumi- 
nosity absorbed and re-radiated by dust in the birth clouds 
and the diffuse ISM. We use the model library described in 
Ida Cunha et all (120081) . which includes 50000 attenuated stel- 
lar spectra spanning a wide range of star formation histories, 
metallicities and dust optical depths. 

3.2. Infrared to millimetre emission 
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In t he context of the model described in Ida Cunha et al.l 
(2008|), the energy absorbed by dust is re-radiated at infrared 
wavelengths through four different components: 

(i) the emission by polycyclic aromatic hydrocarbons 
(PAHs); 

(ii) a hot mid-infrared continuum (with temperature in the 
range 130 -250 K); 

(iii) emission by warm dust in thermal equilibrium (with 
temperature in the range 30 - 60 K and dust emissivity 
index (i = 1.5); 

(iv) emission by cold dust in thermal equilibrium (with tem- 
perature in the range 15 - 25 K and dust emissivity index 
P = 2). 

Ida Cunha et "all (120081) use a wide library of infrared emission 
spectra where the temperatures and relative contributions of 
each component to the total infrared luminosity span a wide 
range of values. The models in this library are then directly 
compared to infrared observations of galaxies to constrain 
each dust emission component. In this paper, since we hardly 
have any observational constraints on the infrared SED of the 
galaxies, we do not require such a wide range of models. 
Therefore, we will adopt a reduced set of dust emission mod- 
els which reflect the range of infrared SED shapes of local, 
normal star-forming galaxies. 

For this work, our goal is to obtain a range of pos- 
sible infrared SEDs that are consistent with the observed 
optical/near-infrared emission in terms of their overall energy 
balance. Therefore, for simplicity, we fix the values of most 
free parameters controlling the shape of the infrared SEDs of 
galaxies t o those of three repres entative infrared SEDs pre- 
sented in Ida Cunha et al.l (120081) : a 'standard' model (with 
equilibrium temperatures of the cold and warm dust compo- 
nents 22 and 48 K, respectively), a 'hot' model (25 and 55 K) 
and a 'cold' model (18 and 40 K); the relative contributions 
of the cold and warm dust components, as well as the PAHs 
and hot mid-infrared continuu m are different for eac h model 
and are listed in Table 1 of Ida Cunha et al.l (120081) . These 
three models were calibrated using observed IRAS and ISO 
infrared fluxes of local star-forming galaxies and span the 
range of observed infrared colors for this low-redshift sam- 
ple. The 'standard' model reproduces the median infrared 
colors of local galaxies, the 'hot' model is representative of 
the warmest observed infrared colors, and the 'cold' model is 
representative of the cooler infrared colors. We build a simpli- 
fied dust emission model library in which we fix the values of 
the dust temperatures and the contribution by PAHs, hot mid- 
infrared continuum and warm dust to the total luminosity of 
birth clouds, as well as the contribution of cold dust to the total 
dust luminosity of the diffuse ISM, to the values of these rep- 
resentative models, while leaving the fraction of total dust lu- 
minosity contributed by the diffuse ISM, = I^St /-^dust. 
and the total dust luminosity, Ldust as free parameters. 

3.3. Radio continuum 

In addition to thermal dust emission, the (sub-)millimeter 
continuum emission of star-forming galaxies can have a non- 
negligible radio continuum emission component, especially 
at the lowest frequencies. This emission is mainly free- 
free emission from H II regions and synchrotron radiation 
from relativistic electrons accelerated in supernova remnants 



(e.g. lCondonlll992l) . Since the Ida Cunha et all d2008l) models 
do not include radio emission, we add a radio continuum com- 
ponent to our model SEDs in order to account for this extra 
contribution to the (sub-)milli meter continuum. We use the 
simple prescription described in D ale & Heloul (120021) . which 
is based on the observed radio/far- infrared correlation. The 
radio/far-infrared correlation (e.g. iHelou et aTlll985l ICondonl 
Il992t lBellll2003l) implies that the ratio between the observed 
far-infrared flux of a galaxy (between 42.5 /im and 122.5 ^tm) 
and the radio flux density at 1.4 GHz, q, is constant. We model 
the radio continuum in star-forming galaxies as a sum of two 
main components: 

(i) a thermal component, consisting mainly of free-free 
emission from ionized gas, with spectral shape / th oc 

V -Q.l. 

(ii) a non-thermal component, consisting of synchrotron ra- 
diation, with spectral shape /„ nth oc i/~ 08 . 

In normal star-forming galaxies, the contribution of the ther- 
mal (free-fr ee) component to the radio continuum at 20 cm 
is ~ 10% (Condon 1992). This allows us to fix the spec- 
tral shape of our radio continuum, which we normalize rela- 
tive to the far-infrare d flux of each m odel in our library using 
the value found by lYun et al] (120011) , q = 2.34. This pre- 
scription is based on the assumption that the galaxies fall in 
the observed radio/far-infrared correlation, and has the advan- 
tage of not requiring any extra free parameters to estimate the 
radio continuum. We also assume that the radi o/far-infrared 
correl ation remains constant with redshift (e.g. Sarge nt et al.l 

HoToh . 

3.4. SED fitting method 

The libraries of attenuated stellar emission and dust emis- 
sion are combined by associating models with similar values 
of — X^ust /-^dust (the fraction of total dust luminosity 
contributed by the diffuse ISM), which are scaled to the same 
total dust luminosity Ldust- This ensures, for each model, the 
energy balance between the radiation absorbed and re-emitted 
by dust in the diffuse ISM and stellar birth clouds. 

For each galaxy, we compare the observed optical fluxes 
in the ACS and NICMOS bands (Section |2| to the predicted 
fluxes for every model of the stochastic library described 
above, by computing the \ 2 goodness of fit for each model. 
We then build the likelihood distribution of any given physical 
parameter for the observed galaxy by weighting the value of 
that parameter in each model by the probability cxp(— x 2 /2). 
Our final estimate of the parameter is the median of the like- 
lihood distribution, with an associated confidence interval 
which is the 16th-84th percentile range of the distribution 
(this confidence interval is tighter for well-constrained param- 
eters and wider when the parameters are not well constrained 
by the available observations). In what follows, the values of 
the physical properties of the galaxies mentioned refer to the 
median values of the probability density distribution. We also 
obtained the best-fit model SED for each galaxy, which is the 
model that minimizes \ 2 ■ 

4. PREDICTED (SUB-)MILLIMETER PROPERTIES 

4.1. SED fitting outputs 

We use the method described above to fit the observed pho- 
tometry of each galaxy and produce likelihood distributions 
of their stellar mass, star formation rate, dust optical depths, 
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Figure 2. Example of a fit to the observed optical spectral energy distribution of a UDF galaxy at z = 1.867 (red points). The black line represents the best-fit 
model SED; the blue line represents the unattenuated (i.e. dust-free) stellar spectrum. The grey-shaded area represents the range of all infrared dust emission 
models in our model library that are consistent with the observed optical/near-infrared fluxes in terms of energy balance. The orange points represent the median 
of the probability density function (PDF) for the predicted fluxes at 345, 230 and 100 GHz (three arbitrary chosen (sub-)mm bands), and the associated error bars 
represent the confidence range, i.e. the 16th-84th percentile range of the PDF. The residuals of the fit are plotted in the panel at the bottom of the SED. The 8 
bottom panels show the PDFs of several parameters: star formation rate; stellar mass; total dust luminosity; dust mass; far-infrared luminosity between 42.5 and 
122.5 (im; and flux densities in the three randomly chosen bands at 345, 230 and 100 GHz. 



infrared luminosities, dust masses, continuum and molecular 
line fluxes in the (sub-)millimeter range. 

In Fig. |2] we show an example SED fit and the associ- 
ated likelihood distributions of some physical parameters: the 
star formation rate (SFR), stellar mass (M»), dust luminos- 
ity (LdustX dust mass (AidustX far-infrared luminosity (£fir, 
defined as the integral of the infrared emission between 42.5 
and 122.5 fim), and the continuum flux densities in a num- 
ber of accessible (sub-)millimeter windows at 345, 230 and 
100 GHz (specifically, ALMA bands 7, 6 and 3, respectively, 
and PdBI/NOEMA band 1, 3 and 4). The ultraviolet to near- 
infrared part of the SED is the stellar population model that 
best fits the data. The far-infrared and (sub-)millimeter part 
of the SED that is plotted corresponds to the model with the 
best fit idust an d /u> but all the other parameters control- 
ling the shape of the SED at these wavelengths are randomly 
drawn from the library of dust emission models. The grey 
shaded area shows the range of possible infrared SEDs al- 
lowed within the uncertainties in dust luminosity with differ- 
ent dust temperatures and contributions of PAHs, mid-IR con- 



tinuum and dust in thermal equilibrium reflecting the diversity 
of possible dust emission models in our model library. The 
orange points show the (exemplary) median-likelihood esti- 
mates of the flux densities at 345, 230 and 100 GHz, where the 
error bars (16th - 84th percentile range) reflect all the differ- 
ent combinations of infrared models that are consistent with 
the observed optical data. 

We impose a minimum of three photometric bands to define 
the SED, and we discard galaxies with z > 5 and galaxies for 
which the fit residuals are larger than 2<r (where a is the un- 
certainty associated with the flux measurement) in each band. 
These criteria allow us to discard the most unreliable SED 
fits: at very high redshift, our model becomes uncertain and 
the observations sample only the far-UV emission, making it 
very difficult to constrain the SED; also, galaxies with very 
high residuals in a given band may indicate problems with the 
photometry or wrong photometric redshift, or the presence 
of an AGN (which is not included in our models). Our final 
sample used in the remainder of this paper consists of 13, 099 
sources with z < 5. 
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4.2. Derived physical properties 

Stellar masses are constrained to within ±0.35 dex, which 
reflects uncertainties due to the fact that, for most of our 
galaxies, observations do not include the rest-frame near- 
infrared, where the light is dominated by low-mass stars, 
which constitute the bulk of the stellar mas s in g alaxies. How- 
ever, we show in the next section (Section 14. 31 > that this does 
not cause any systematic effects on the stellar mass estimates. 
The star formation rates are constrained to within typically 
±0.2 dex, due to the fact that the observed SEDs sample the 
emission by young stars in the rest-frame ultraviolet. The 
dust luminosities are more uncertain (confidence ranges are 
typically ±0.45 dex), as expected due to the lack of infrared 
observations for our sample. The dust luminosity is estimated 
by our model by calculating the total energy absorbed by dust, 
taking into account the light emitted by stars and the attenua- 
tion by dust. Therefore, by construction, our dust luminosities 
are consistent with all stellar and dust attenuation parameters 
(SFR, M*, [i, fy, /„,), from an en ergy balance perspective 
(as described in da Cunha et al. 2008 and in Section|3]l. Even 
though we have significant uncertainties in the dust luminos- 
ity estimates, as expected from our sparse SED sampling, we 
are still able to predict Ldust well within an order of magni- 
tude (we analyze possible systematic effects in Section |4~3"T >. 
The dust masses are also estimated by using all the dust emis- 
sion model templates that are consistent with the statistical es- 
timates on dust luminosities. The confidence range for Md us t 
is very large (±0.55 dex), and reflects not only the uncertainty 
in £dust, but also the large uncertainty in infrared SED shapes 
and dust temperatures. These dust masses are merely indica- 
tive of the range of dust masses that are consistent with the 
observed SEDs in terms of energy balance, and taking into 
account a range of possible dust temperatures, 18-25 K for 
the cold dust (with f3 — 2), and 40-55 K for the warm dust 
(with/3 = 1.5). 

The distributions of physical parameters inferred from our 
SED fits (Fig. show that the bulk of our galaxies are low- 
mass, low star formation rate and low dust attenuation (fy 
and jitfy) sources, i.e. blue star-forming galaxies (consistent 
with the finding of a l arge population of faint blue galaxies in 
the UDF described in lCoe et afll2006l) . As expected, galaxies 
in the highest redshift bin, z > 2.5, have typically higher stel- 
lar masses and star formation rates, because only the bright- 
est galaxies are detected. The median dust luminosity of the 
sources, Ldust increases from log(Xdust/£o) — 8.0 in the 
lowest redshift bins to log(idust/i©) — 9.2 at z > 2.5. 
Figsj3jg) and Oh) show that this is not necessarily due to an 
increase in the dust optical depth of galaxies in the highest 
redshift bin, but rather to the fact that the galaxies detected 
have larger stellar masses and star formation rates, as shown 
in Figs|fla) and[2b). 

4.3. Reliability of infrared luminosity estimates 

It is important to test how well we can predict the total 
infrared luminosity of the galaxies in our sample from their 
observed rest-frame UV/optical spectral energy distributions. 
To do so, we use the sample of 54 UDF galaxies detected 
in the Herschel/PACS and Herschel/SPIRE bands as part of 
the GOODS-Herschel program described in Section |2~2l For 
each of these galaxies, we have well-sampled SEDs from the 
ultraviolet to the far-infrared, which allows us to test our SED 
extrapolations. We fit the observed (more complete) SEDs 
of this subsample of galaxies using the same method as de- 



scribed in Section l3~4l First, we include the full observed ul- 
traviolet to far-infrared observations in the SED fits, not only 
to check that our model can reproduce consistently the SEDs 
of these galaxies, but also to get the best possible estimates 
of the stellar masses, star formation rates, dust luminosities 
and continuum (sub-)mm fluxes for these sources. Then, we 
re-fit the SEDs using only observations between the U band 
and the K band, to mimic the set of observations available for 
the majority of galaxies discussed above. 

In Fig.@Ja), we compare the median-likelihood of the Her- 
schel/SPIRE 250-/im flux derived from the fits from the U 
to K band with the actual observed 250-/im flux for each 
galaxy. We find a small systematic offset of 0.25 dex between 
the observations and our estimates, in the sense that we tend 
to underestimate the 250-/im flux of the galaxies on average 
with our U-lo-K-bwd fits. However, for most of the galax- 
ies, the observed value is still within the confidence range de- 
rived from our fits. This effect is likely to be less important 
for the total dust luminosity and (sub-)mm fluxes, as these do 
not depend as strongly on the exact location of the peak of 
the infrared SED (i.e. the actual dust temperature spanned by 
the models). In the next three panels of Fig. H] we compare 
the median-likelihoods of the stellar masses, dust luminosi- 
ties, and continuum ALMA Band 6 fluxes obtained with the 
two sets of SED fits (CZ-to-K-band fit vs. [/-to-SPIRE-fit). 
Fig-Etb) shows that the stellar masses agree remarkably well 
between the two fits (with a dispersion around the identity line 
of 0.13 dex). Not surprisingly, the total dust luminosities and 
predicted ALMA band 6 fluxes do not agree as well, as shown 
in Figs. Uc) and (d). The inclusion of infrared data in the fits 
helps constrain these properties better, as shown by the signif- 
icantly reduced confidence ranges. In the case of £dust> this 
happens because the infrared data allow us to constrain the 
bolometric dust luminosity by fitting the dust emission itself 
(as opposed to constraining idust from the attenuated spec- 
trum alone); the constraints on SFR are also tightened because 
we can account for dust-obscured star formation rate more ac- 
curately. The different far-infrared fluxes obtained with PACS 
and SPIRE help constrain the shape of the dust SED, namely 
the dust temperatures and relative contributions of the warm 
and cold dust components to the total infrared emission. This 
helps obtaining tighter constraints on the (sub-)mm contin- 
uum fluxes (namely the ALMA Band 6 flux shown as an ex- 
ample in Fig.|4fd)). 

Even though the inclusion of infrared data helps constrain- 
ing parameters such as the star formation rate, dust luminos- 
ity and ALMA continuum fluxes, the median likelihood es- 
timates of these parameters when excluding the infrared data 
agree very well with the estimates derived from the full SED 
fits, even if, as expected, the associated confidence ranges 
are larger. We find very small offsets between the averages 
of the median-likelihood estimates derived from the two fits: 
0.02 dex for SFR, 0.08 dex for L dus t, and 0.07 dex for the pre- 
dicted ALMA Band 6 continuum flux (in the sense of the U- 
to-fC-band fit slightly underestimating the parameters), with 
a dispersion of ~ 0.40 dex for all cases. These very small 
systematic offsets are well within our fit confidence ranges, 
and show that our approach to predict infrared luminosities 
and (sub-)mm continuum fluxes from modelling UV/optical 
SEDs is reliable. We note however that the difference be- 
tween 250/im fluxes and total dust luminosities derived from 
the fits to the UV/optical data only and those measured with 
Herschel correlates with the dust optical depth in the galax- 
ies. We tend to underestimate the (sub-)mm fluxes/dust lu- 
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Figure 3. Distributions of the median-likelihood parameters derived from SED fitting for the whole Hubble UDF galaxy sample, (a) stellar mass; (b) star 
formation rate averaged over the last 100 Myr; (c) specific star formation rate, defined as the star formation rate divided by stellar mass; (d) total dust luminosity; 
(e) dust mass; (f) fraction of total dust luminosity contributed by the diffuse ISM, f^; (g) effective V-band optical depth seen by stars younger than 10 Myr in 
birth clouds, fy ; (h) effective V-band optical depth seen by stars older than 10 Myr in the diffuse ISM. The grey histograms represent the whole sample, and the 
colored histograms represent the distribution of parameters of galaxies divided in three redshift bins: green: z < 1.5; yellow: 1.5 < z < 2.5; red: z > 2.5. The 
error bars in the top right-hand corner of each plot represent the median confidence range for each parameter. We note that the sharp drop towards lower values 
of stellar mass, SFR and dust mass/luminosity is due to incompleteness of the photometric catalog towards fainter flux levels (see Fig.[TJ. 



minosity when using only the [7-to-i\~-band fits for galaxies 
with the highest dust attenuations (which translate into high 
infrared-to-optical ratios). Thi s is due to the fact th at our dust 
attenuation prior (Section [3] Ida Cunha et al.l 12008b leads to 
an underestimation of the optical depth for extremely dust- 
enshrouded, starburst-lik e sources (such as loc al ULIRGs or 
high-redshift SMGs; see Ida Cunha et ail 120101) . While these 
galaxies can be a negligible fraction of our sampl e in num- 
ber (e.g. lRodighiero et al.ll201 UlSargent et al.ll2012|) , they can 
dominate the bright (sub-)mm counts. This is the case for the 
two GOODS/Herschel sources in our sample with the highest 
redshift, which are marked in Fig. |4]with squares. Due to the 
flux limit of this sample, at the highest redshifts (z ~ 2), only 
very dust-obscured ULIRG-type galaxies were selected. For 
this particular type of galaxies, the SED models used in Sec- 
tion[3]become limited. However, we expect this kind of galax- 
ies to be rare in our optically-selected sample of the UDF, 
and therefore we do not expect this limitation to greatly affect 
our results. We also note that our optically-selected catalogue 
is also likely to mis s completely optic ally-obsc ured galax- 
ies (e.g. HDF850.1. IWalter et al.ll20Tl GN10. iD^ddieFaT] 



I2009at GN20. lDaddi et alj|2009bl) . The very good agreement 
between parameters derived from fits to the UV/optical SED 
versus parameters derived from fits to the full SED is con- 
sistent with previous results showing that the star formation 
properties of normal star-forming galaxies up to z ~ 2 can 
be reliably derived from UV/optical obser vations alone (e.g. 
iDaddi et alf2005l 120071: liteddv et alj|2006l) . This implies that 
the ISM of these galaxies is not heavily optically-thick, and 
we can apply our energy balance technique to interpret the 
SEDs of most normal star-forming 'main sequence' galaxies. 

4.4. Number counts 

As a consistency check, we now compare our continuum 
flux density predictions with previously obtained number 
counts at (sub-)millimeter wavelengths. 

Number counts at 850/im have been obtained using 
the SCUBA bolometer array on the JCMT and LABOCA 
on APEX by a number of groups ov e r the last decade 
(e.g. IScott etal] 12001 ISmail et alj|200l iBorvs et ail [2003; 
Coppi n et al.l 120061: IScott et al.l l2006t iKnudsen eTail 12008 ; 
WeiB et aLll2009h . These counts can be directly compared to 
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Figure 4. (a) Comparison between the observed Herschel/SPJRE 250-/tm flux (x-axis) and our Bayesian median-likelihood estimate of the 250-fim flux of 
each galaxy based on SED fits from the £/-band to the X-band (y-axis). The other three panels show the comparison between our Bayesian median-likelihood 
estimates of GOODS/Herschel galaxy parameters obtained when fitting the full SED from the U-band to the longest available SPIRE band (x-axis), and when 
fitting only the SED from the {/-band to the K-band (y-axis): (b) stellar mass; (c) total dust luminosity; (d) continuum flux in the ALMA band 6 (230 GHz). 
Each galaxy is color-coded according to redshift. The error bars show the 16th - 84th percentile range of the likelihood distributions. In all panels, the grey solid 
line is the identity line, and the dotted lines show a ±0.5 dex offset for reference. The two points marked with crosses are galaxies that show a significant AGN 
contribution in the infrared; the two points marked with squares are galaxies which show a ULIRG-like SED, i.e. they seem to be very optically thick (given their 
high intrinsic infrared-to-optical emission ratios). Our SED modelling may not be reliable for these four galaxies, but overall we find a good agreement between 
the estimates derived from fitting the full SED and those from fitting the SED only up to the iC-band. 



our predictions at 345 GHz (ALMA band 7, PdBI/NOEMA 
band 4) since this band probes roughly the same wavelength. 
In Fig. [5] we compare our predicted cumulative number 
counts at 345 GHz with these previous observations. Our pre- 
dicted number count range to first order agrees with the ob- 
served number counts, but we do not reach the higher fluxes 
probed by sub-mm observations. The lack of the brightest 
sources is due to two reasons. First, the field on which we 
base our predictions is very small (the size of the UDF is 
only 3.3 x 1CP 3 deg 2 ), i.e. we do not expect the presence 
of a significant population of bright sources in the field. Sec- 
ond, we are working with an optically-detected sample, and 
the bright sub-mm counts are dominated by optically thick 
sources which are likely not detected in the optical. For ex- 
ample, LESSJ0333243. 6-274644, the only sub-mm source 
detected in the Hubble UDF as part of the LESS survey 
(LABOCA observations of the Extended Chandra Deep Field 
South ), with a flux density of 6.4 mJy at 870/^m dWeifi et alj 
2009), has no optical counterpart in our catalog, presum- 
ably because it is an optically-thick sub-mm galaxy (SMG; 



cf. lDunlopir201 11) . We find that the number counts at fluxes 
> 0.5 mJy are dominated by galaxies with fify > 1, i.e. 
where the ISM is optically-thick on average. 

As an additional test on our number count predictions, we 
turn to a wider field covered by the LESS survey. To do so, we 
expand our analysi s to the FIREWORKS photometric catalog 
dWuvts et alj|2008l) on the CDF-S, described in Section |2T2l 
FIREWORKS covers an area that is about 10 times the area 
of the UDF and about 10 times smaller than the full E-CDFS. 
The photometric catalog is much shallower than that for the 
UDF. For the area covered by FIREWORKS, we estimate 
between 6 and 23 sources to have 870-/im fluxes above 4.7 
mJy (the flux li mit of the LESS catalog). For comparison, 
iWeiB et all (120091) find 10 sources over the same area. Our 
prediction is thus broadly consistent with the LESS measure- 
ments. 

4.5. Extragalactic background light 

We can also compare our predictions with measurements 
of the integrated extragalactic background light in the sub- 
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Figure 5. Predicted cumulative number counts in Band 7 at 345 GHz (black solid line), with confidence range estimated using the upper and lower flux limits 
given by the confidence range for each galaxy (dotted lines). The colored point s show observ ed values at ~850/im from different studies (see figure legend for 
references). The vertical line shows the flux limit of the LESS survey catalog ( WeiB et al. 2009), which includes the UDF. 



mm. Using our median-likelihood estimate of the 345 GHz 
flux density for each galaxy in the Hubble UDF, we ob- 
tain an integrated continuum 870-/im flux density of 45.6 Jy 
deg" 2 . If we add the contribution from LESSJ0333243.6- 
274644 (which is not part of our sample but is detected in the 
LESS survey with a 6.4 mjy flux), we obtain an EBL value 
of 47.5 Jy deg~ 2 . This value is fully consistent with mea- 
surements of the extragalactic background light from COBE 
obser vations 44± 15 Jy deg~ 2 (Pu get et al.119 96; Fixse n et alj 
[19981) . We note that the use of fixed spectral energy distribu- 
tion templates to derive the (sub-)millimeter continuum fluxes 
of the galaxies would lead to EBL values that are inconsistent 
with the observed value (see Appendix lAt. In TableQ] we list 
our estimates of the EBL in different (sub-)millimeter bands. 

The estimated EBL at 870 fjm using the FIREWORKS cat- 
alog over the CDF-S area is 36.3 Jy deg -2 , broadly consistent 
with our estimate based on the UDF area only (Table[T). This 
is lower than the COBE observed value quoted above, but it is 
still consistent with the observations, since the FIREWORKS 
catalog does not reach very deep, and therefore it is likely to 
miss the large number of faint galaxies that make up for a 
significant fraction of the extragalactic background light. 

4.6. Continuum flux density predictions 

In this section, we present our general predictions for the 
continuum flux densities of our galaxies. In Table[6] we make 
our continuum predictions in all relevant (sub-)mm bands 
available for all the galaxies in our sample. 

As an example, in Fig. [6] we plot continuum map of the 
UDF at 230 GHz (ALMA band 6, PdBI/NOEMA band 3) 
using our predictions (right-hand panels). Such images can 
be directly compared to future deep fields performed with 
ALMA or other facilities. We note that this figure is based 
solely on our optically-based predictions, and so they are 
missing the only known bright (sub-)mm source in the UDF: 
LESSJ0333243.6-274644. The top panels of Fig.[6]show our 



Table 1 

Predicted extragalactic background light. 



Frequency 


Observatory 


EBL 


/GHz 




/ Jy deg -2 


38 


ALMA 1, JVLA Ka 


0.10 


80 


ALMA 2 


0.74 


100 


ALMA 3, PdBI/NOEMA 1 


1.53 


144 


ALMA 4, PdBI/NOEMA 2 


4.48 


230 


ALMA 6, PdBI/NOEMA 3 


18.2 


345 


ALMA 7, PdBI/NOEMA 4 


45.6 


430 


ALMA 8 


63.1 


660 


ALMA 9 


123 


870 


ALMA 10 


146 


Note. — 


Estimates of the extragalactic background 



light at different frequencies using our flux predictions 
(area of the UDF field is 0.0033 deg 2 ) for the 13,099 
galaxies in our sample. 



full sample, and we then divide the sample in two redshift bins 
( z < 1; middle panels) and (z > 1; bottom panels). This illus- 
trates the differences between galaxy detections as a function 
of redshift between the optical and the (sub-)mm, in partic- 
ular that we expect galaxies to be relatively brighter in the 
(sub-)mm at high redshift thanks to the negative k-correction. 
Therefore, we will be able to detect 'normal' galaxies out to 
higher redshifts in the (sub-)mm with new facilities, as we 
discuss in more detail in Section 

In Fig. [JJ we plot the distribution of the predicted contin- 
uum flux densities of all the galaxies in our sample in all 
current and future ALMA, JVLA and IRAM PdBI/NOEMA 
bands from 38 GHz to 870 GHz. We plot the expected num- 
ber of galaxies per flux bin in the total UDF area in each band. 
The distribution of fluxes peaks at higher fluxes in the highest 
frequency ALMA bands, because, even taking into account 
fc-correction effects, these bands sample the emission from 
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HST/ACS V-band (whole sample) 



continuum at 230 GHz (whole sample) 



45' 



46' 



47' 



48' 



49' 



-27°50' 



45' 



46' 



47' 



48' 



49' 



-27°50' 



.', „• ••. 

' ' A-.. V * 'I 



03 h 32 m 48 s 42 s 39 s 36 s 33 s 30 s 
HST/ACS V-band (z < 1) 



• f ... ■ . . 

* • « ' . - 



03 h 32 m 48 s 42 s 39 s 36 s 33 s 30 s 



HST/ACS V-band (z >= 1) 



-27°50 




03 h 32 m 48 s 42 s 39 s 36 s 33 s 30 s 



45' 



46' 



47' 



48' 



49' 



■27°50' 



45' 



46' 



47' 



48' 



49' 



■27°50' 



45' 



46' 



47' 



48' 



49' 



27°50' 



-:;r.:-.. 



03 h 32 m 48 3 42 s 39 s 36 s 33 3 30 s 
continuum at 230 GHz (z < 1) 



03 h 32 m 48 3 42 s 39 s 36 s 33 3 30 s 
continuum at 230 GHz (z >= 1) 



•-. ».*. 

.. »•<; , 
v -a... :*: *»•" 

'.•>•*•. 

* . . •.•••.%,.« - 

■• ••»•'•>....••■' : ' 5 

i . •• \ /-V. • ■-*•■> ". 
tv . r. : * 
•. > • !?»• .v. 



03 h 32 m 48 s 



42 s 39 s 36 s 33 s 30 s 









c 


0.8 


\ 








< 


0.6 


X 




o 






0.4 




0.2 





Figure 6. Comparison between the UDF observed in the V-band (left-hand panels) and at 230 GHz (right-hand panels; using our predictions); top panels: whole 
sample; middle panels: galaxies with z < 1; bottom panels: galaxies with z > 1. The V-band image was generated using ACS F606W image of the UDF 
(from the HST archive); the 230-GHz image was generated assuming point sources and convolving with the ALMA synthetic beam in the compact configuration, 
1.5 arcsec. The greyscale bar shows the predicted (sub-)mm fluxes of the galaxies; no noise is included in either the optical or (sub-)mm maps. We note that 
these maps do not include the SMG galaxy detected in the UDF as part of the LESS survey (LESSJ0333243.6-274644). 
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Figure 7. Distribution of the predicted continuum flux densities per flux density bin in different (sub-)millimeter bands for all the galaxies in our sample. The 
frequency is indicated in the bottom-left corner of each panel, as well as the corresponding bands in different observatories. The median confidence range for the 
continuum flux in each band is plotted in the upper right-hand corner of each panel. 



12 



E. DA CUNHA ET AL. 



the galaxies closer to the peak of the dust SED. In Section[3] 
we discuss the feasibility of performing a blind survey of the 
UDF with ALMA at full operation, and use these predicted 
fluxes, combined with the projected ALMA sensitivities, to 
obtain an estimate of the expected number of continuum de- 
tections. 

4.7. CO and [CII] line predictions 

ALMA and JVLA will detect CO and [CII] line emis- 
sion from high-redshift galaxies, which will allow us to 
determine redsh ifts, molecular gas reservoirs, and dynami- 
cal masses ("e.g. [Solomon & Vanden Bo ut 2005; Dad di et alJ 
20l0i IGenzel et all 120 1 Ot ITacconi et all BO 1 Ot IWalter et all 
201 lb . The rest-frame frequency, v res u of the CO lines corre- 
sponding to J — > J — 1 transitions from J = 1 to J = 7 and 
of the [CII] line are given in Table[5] The observed frequency 
of each line varies with redshift as t'obs = ^ rC st(l + in 
Fig. [T4l (Appendix IB1. we plot the observed frequency of the 
seven first CO transitions, CO(l-O) (i.e. J = 1) to CO(7-6) 
(i.e. J = 7) and of the [CII] line as a function of redshift, with 
the frequency ranges covered by each ALMA, PdBI/NOEMA 
and JVLA band shaded in grey and green. In Table [5] (Ap- 
pendix[JU, we explicitly list the redshift ranges where the CO 
lines and [CII] are observable in each band. It is clear that 
the lowest frequency bands, such as JVLA K, Ka and Q will 
be crucial to probe the low J CO transitions, in particular the 
CO(l-O) line, in z > 1 galaxies. All the other PdBI/NOEMA 
and ALMA bands will potentially detect higher J CO transi- 
tions at different redshifts, depending on the excitation state 
of the gas in galaxies. The highest-frequency ALMA bands 
will not only sample high- J CO lines at low redshifts and the 
continuum dust emission nearest to its peak, as mentioned in 
the previous section, but also the [CII] line out to high red- 
shifts. 

In this section, we attempt to predict the CO and [CII] line 
fluxes for the galaxies in our UDF sample using empirical re- 
lations that relate line luminosities with the infrared luminos- 
ity of the galaxies, for w hich we have a statistical estimate 
from our SED fits (SectionHTTTi. 

4.7.1. CO emission 

The CO line luminosity of galaxies depends on vari- 
ous factors, such as the gas heating by starbursts, AGN, 
and the cosmic microwave background at high redshifts 
(e.g. lCombes et al.lll999trObreschkow et alj|2009ri da Cunha 
et al., in pre p.), as well as the clumpin ess and metallicity of 
the gas (e.g. Obreschkow et al. 2009b). In this section, for 
simplicity, we predict the CO line luminosity of the galaxies 
in our sample using simple, empirically calibrated prescrip- 
tions. It has been found for a wide range of galaxy types, both 
in the local and high-redshift Universe, that the CO line lumi- 
nosity of star-fo rming galaxies correlates with their infrared 
luminosity (e.g. Solomon & Vanden Boutll2005l IGenzel et alJ 
l2010l IDaddi et alj|2010bh ~ 

The following relation between C O line luminosity and far- 
infrared luminosity was derived by iDaddi et al.l (1201 Obi) for 
BzK galaxies, i.e. gas-rich star-forming disks at high red- 
shifts: 

log(i m ) = 1.13 logOL'co) + 0.53 , (1) 

where Lir is the infrared luminosity (integrated between 8 
and 1000 /im) in Lq and L' co is the CO(1-0) line lumi- 
nosity in K km s _1 pc 2 . We obtain an estimate of L' co 



using this equation and the statistical estimate on Lm ob- 
tained for the SED fits of our galaxies; the resulting distri- 
bution of L' co for the whole sample is plotted in the top left- 
hand panel of Fig. [8] We chose this empirical calibration be- 
tween Lin and L' rn because our physical parameter estimates 
in Section |4~2l indicate that most of these galaxies would be 
comparable to normal, 'main-sequence' star-forming disks, 
with typical infrared luminosities Lm < 10 11 Lq. We note 
that eq . Q] is similar to the relation found by Genz eTet al.l 
(2010) for isolated, star-forming galaxies out to z = 2. 
Other calibrations of this relation have been derived which 
include more extreme galaxies such as starbursts, mergers, 
and A GN (e.g. lSolomon & Vanden Boutll2005l;lRiechers et all 
2006). When including these extreme galaxies, the relation 
betwe en infrared and CO line l uminos ity becomes steeper, 
e.g. Solomon & Vanden Bout (2005) find log(LFiR.) = 
1.71og(i^ ) — 5.0. Using this steeper relation for the in- 
frared luminosity range of our galaxies (Ldust < 10 11 Lq) 
would result in CO line luminosities over one order of magni- 
tude higher than those predicted using Eq. [T]for the Lm range 
of the galaxies in our sample. We discuss the implications of 
using these different assumptions for the predicted number of 
CO line detections in Section[5] 

From L' co (computed using Eq. [TJ, we then compute the 
corr esponding flux of the CO(l-O ) line, S„ , using 

(e.g. lSolomon & Vanden Boutll2005l) : 

L' co = 3-25 x 10 7 ^o(i-o) Av „-* ^ (1 + z) -3 ) (2) 

where 5?°^ °' ) is the flux density in Jy, Av is the line 
width in km s _1 (the velocity-integrated flux of the line is 
ico(i-o) = S v Av), ^ bs is the observed frequency 

of the line in GHz, and Dl is the luminosity distance in 
Mpc. We assume a typical line width of 300 km s _1 , consis- 
tent with typical line- widths measured in high-redshift star- 
forming galaxies (e .g. IDaddi et"alll2010bt IGenzel et al.ll2010T: 
ITacconi et al.M2.Q10l) . In the top right-hand panel of Fig. [8] 
we plot the distribution of the velocity-integrated flux of the 
CO(1-0) line for all the galaxies in our sample computed us- 
ing equation [2] 

The fluxes of higher transition CO lines depend highly on 
the excitation of the CO gas in galaxies. Different physi- 
cal conditions in the gas produce different CO spec tral line 
energy distributions (SLEDs; e.g. Wei ss et all 120071) . which 
translate into different ratios between the CO(1-0) line and 
the higher J lines. To compute the predicted flux of the 
CO(2-l), CO(3-2), CO(4-3), CO(5^), CO(6-5) and CO(7- 
6) line s, we assume two possible CO SLEDs from Weis s et all 
(2007): the Milky Way CO SLED and the M 82 center CO 
SLED. These two cases correspond to very low and high ex- 
citation of the gas, respectively, and should bracket a large 
realistic range of possible physical conditions in star-forming 
galaxies. In the six bottom panels of Fig. [8] we show the dis- 
tribution of expected velocity-integrated CO line fluxes for the 
galaxies in our sample, and compare the predictions for these 
two excitation scenarios. For each CO line, these two extreme 
excitations should bracket the full range of line fluxes ex- 
pected for our sample of star-forming galaxies (as supported 
by observations of multiple CO lines in a wide range of sys- 
tems fromJ_o£aLquiescent galaxies to high-redshift QSOs; see 
e.g. iWeiss et al.M2.Q07h . The difference between the predic- 
tions of CO line fluxes using these two CO SLEDs increases 
with increasing J: the higher- J CO lines are stronger in the 
case of the M 82 center SLED, which corresponds to a higher 
CO excitation. In Table |7l we provide the predicted CO line 
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Figure 8. Top panels: distribution of total CO line luminosity L' C q (top left) and velocity-integrated flux of the CO(1-0) line (top right) for all our sample. Other 
panels (stalling at second row): distribution of the predicted velocity-i ntegra ted CO line fluxes for transitions J = 2, 3, 4, 5, 6, and 7. These fluxes are computed 
from the predicted CO(1-0) line, using the CO SLEDs of Weiss et al. ( 2007) for the Milky Way (light gray histogram) and the center of M 82 (dark gray, dashed 
histogram). These two extreme CO SLEDs should bracket a large realistic range of possible CO excitations of star-forming galaxies (e.g. W eiss et 
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fluxes for all the galaxies in our sample (including both CO 
excitation scenarios). 

To predict the number of expected CO line detections given 
a certain flux limit in various (sub-)millimeter bands (ALMA, 
JVLA or PdBI/NOEMA), we can build the expected distribu- 
tion of CO line fluxes in each band based on the predictions 
described above. First, for each frequency band listed in Ta- 
ble[5] we retain galaxies for which the redshift falls in a range 
where one of the CO lines can be observed (these ranges are 
listed in Table [5j see also Fig. PHI . Then, we compute the 
expected line fluxes, S^ nc in Jy, for each galaxy, by assum- 
ing a typical line width of Av = 300 km s _1 (using eq. |2] 
to the get flux of the CO(1-0) line and the CO SLEDs to get 
the higher- J lines). In Fig.|£l as examples, we plot the dis- 
tribution of line fluxes from CO(1-0) to CO(7-6) if one were 
to fully cover the bands from 201 to 275 GHz (ALMA band 
6, PdBI/NOEMA band 3), 80 to 116 GHz (ALMA band 3, 
PdBI/NOEMA band 1) and 18 to 50 GHz (ALMA band 1 / 
JVLA bands K, Ka and Q). We show the CO line fluxes cor- 
responding to two gas excitation scenarios (i.e. CO SLEDs): 
Milky Way-type (left-hand panels) and M82-type (right-hand 
panels). For example, in the frequency range 80 to 1 16 GHz, 
we plot the distribution of CO(1-0) line fluxes only of galax- 



ies with redshifts z < 0.44, for which the CO(1-0) line would 
be redshifted to that frequency band (Table [5). Similarly, for 
the distribution of CO(2-l) line fluxes in that band, we in- 
clude only galaxies with 0.99 < z < 1.88, and so on until the 
CO(7-6) line. For clarity, the distributions plotted in Fig.|9]are 
cumulative: the darkest-colored histogram shows the number 
of galaxies per line flux bin for the CO(1-0) line, the next, 
lighter histogram shows the number of galaxies per line flux 
bin for the CO(1-0) and CO(2-l) line, the next histogram 
adds the number of galaxies per line flux bin for the CO(3-2) 
line, etc. That is, the lightest-colored histogram shows the to- 
tal number of galaxies per line flux bin when including all CO 
lines from J = 1 to J = 7. In general, Fig. [9] shows that the 
number of galaxies in the highest flux bins is largest if the CO 
SLED is M82-like, as expected, so the number of detections 
given a certain flux limit gre atly depends on the CO excitation 
(as discussed in Section l5T2b . 

4.7.2. [CII] emission 

The [CII] line at 158 /im is the main cooling line of the 
ISM in galaxies, and it arises mainly from photo-dissociation 
regions - at the interface between the ionized gas and the neu- 
tral and molecular gas - which are typically associated with 
star-forming regions. This line is therefore one of the main 
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Figure 9. Number of expected CO line detections per line flux bin for three frequency ranges and two different molecular gas excitation scenarios. The line 
fluxes are computed assuming a line width of 300 km s — . Top: frequency between 201 and 275 GHz (ALMA band 6 and PdBI/NOEMA band 3); Middle: 
frequency between 80 and 116 GHz (ALMA band 3 and PdBI/NOEMA band 1); Bottom: frequency between 18 and 50 GHz (ALMA band 1 and JVLA bands 
K, Ka and Q). The left-hand and right-hand panels assume a Milky Way and M82 CO spectral line energy distribution, respectively. The lowest, darkest color 
histograms show the distribution of CO(1-0) fluxes, the next (lighter-colored) histogram shows the joint distribution of CO(1-0) and CO(2-l) fluxes, i.e. the 
distribution of CO(2-l) can be inferred from the increment in the histogram relatively to the histogram below, and so on until the lightest-colored histogram, 
which shows the sum of the distributions of all the CO line fluxes from J = 1 to J = 7. 



far-infrared tracers of star formation in galaxies (IStacev et al. 
199U iBoselli et al.l l2002t IStacev et aljfeoiOt Ide Looze et at 
201 lb . Since this is the brightest far-infrared line, it will be 
readily detected in deep observations with ALMA, particu- 
larly using the highest-frequency bands. 

We rely on previous observational studies of the ratio of 
[CII] line to far-infrared luminosity (Xrciii / £fir) of galaxies 
to estimate the [CII] line luminosity for each galaxy in our 
sample. The ratio L[ CII ]/Xfir of normal star-forming galax- 



ies varies typically between 1 and 0.1%, and h as been shown 
to ant i-correlate with dust heating intensity (e.g. Brau heret al.l 
2008). Here, for simplicity, and considering the relatively 
large error bars on our Lfir estimates, we adopt a con- 
stant ratio of log(L[cii]/£piR,) = —2.5, which corresponds 
to the average value for normal star-forming galaxies with 
low far-infrared luminosities Lir < 10 11 L© (i.e. simi- 
lar to the galaxies i n our sample) and average dust heating 
(IBoselli et al.H2002t iBrauher et al.ll2008t iGracia-Carpio et all 
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l20llt ICox et alj|20llh . In Fig. Ha), we plot the distribu- 
tion of the expected velocity-integrated flux of [CII], -frciil = 
Si Av for all the galaxies in our sample, computed using 
eq.El 

L [cirj = 1-04 x lQ- 3 Sf^ Av v Iest D 2 L (1 + z)- 1 , (3) 

where Lrcii] is the [CII] line luminosity in L@, Sjp 1 ^ is the 
line flux in Jy, Av is the velocity dispersion in km s _1 , Di is 
the luminosity distance, and v lcst = 1900.54 GHz. At z < 5, 
the observed frequency of the [CII] line falls in the highest fre- 
quency ALMA bands, namely bands 7, 8, 9 and 10 (Eig.fPfli. 
In Fig. flOl' b). we plot the distribution of expected detections in 
these bands as a function of [CII] line flux (computed assum- 
ing = 300 km s _1 ). In Table|7l we provide the predicted 
[CII] line fluxes for all the galaxies in our sample. 

5. POSSIBLE DEEP FIELD STRATEGIES WITH ALMA 

Based on the results presented in the previous sections, we 
now discuss the feasibility of carrying out a deep field survey 
of the Hubble UDF with ALMA at Full Operations (i.e. using 
50 antennas). In the following, we will assume a total time 
of 500 hours for the full survey, and investigate the possible 
setups of such a survey to maximize the redshift coverage and 
number of galaxy detections. 

One immediate drawback of ALMA as a survey instrument 
is that the primary beam size, which is driven by the size of 
the antennas, is relatively small in all bands: the primary beam 
size is given, to first order, by 20". 3 x 300/ (V/GHz), where v 
is the observing frequency. Therefore, even a relatively small 
area field as the Hubble UDF (3.45' x 3.45', i.e. a total area 
of 1 1.9 arcmin 2 ), will be hard to cover with ALMA, and will 
require significant mosaicking. It is beyond the scope of this 
paper to go into the technical details of how such a mosaic 
should be set up. In the following, for illustration purposes, 
we simply consider a total on-source time of 500 hours for the 
UDF, which we divide in different pointings to cover the area 
of the UDF - we do not include overheads and mosaicking 
details (such as degrading sensitivities inside the beam etc.). 

5.1. Continuum detections 

For each ALMA band, given the area of the primary beam, 
we compute how many effective pointings are needed to cover 
the whole area of the UDF (Table|2]i. The number of pointings 
needed to cover the UDF increases from Band 1 to Band 10, 
from only 3 pointings to over 1000 pointings. In this back- 
of-the-envelope calculation we take the integration time per 
pointing in each band as the total 500 hours divided by the 
required number of pointings, and then use the ALMA Sen- 
sitivity Calculator (ASC£] to compute the continuum sensi- 
tivity that can be reached in each integration time, assuming 
a 8 GHz bandwidth and the default weather conditions. In 
Table [2] we list the sensitivities and the expected number of 
3cr continuum detections in each band for the wh ole field, 
based on our continuum predictions of Section 14.61 The re- 
sulting number of predicted continuum detections is a com- 
bination of the intrinsic flux of each galaxy, the fc-correction 
for each galaxy at each redshift, and the changing integra- 
tion time per pointing due to varying beam size; for refer- 
ence, 1 hour integration time corresponds to a rms of 4.81, 
12.6 and 483 /Jy at 38 GHz (ALMA 1), 230 GHz (ALMA 6) 

1 http://almascience.eso.org/call-for-proposals/sensitivity-calculator 



and 870 GHz (ALMA 10), respectively. Table |2] shows that 
the number of expected continuum detections is maximum 
for band 6 at 230 GHz (601 detections), with a large num- 
ber of predicted detections also in band 3 at 100 GHz (221) 
band 4 at 144 GHz (363) and band 7 at 345 GHz (522). In 
Figs.Q]and[TT| we compare the distributions of the properties 
of the 601 galaxies that would be detected at the 3c-level in 
the ALMA band 6 to those of the original sample. We note 
that the redshift distribution of these sources is relatively flat, 
thanks to the negative ^-correction in the sub-mm. Also not 
surprisingly, these figures show that we expect to detect only 
the most massive, highly star-forming and dusty galaxies in 
our sample. These galaxies are still over an order of mag- 
nitude star-forming and dusty than classic SMGs detected in 
blind 'pre-ALMA' (sub-)mm surveys, which have typically 
star formation rates > 100 M Q yr _1 and dust luminosities 
> 10 12 Lq. The median detected galaxy has a stellar mass 
of 3 x 10 9 Mq, a star formation rate of ~ 5 Mqvt -1 , dust 
luminosity of 4 x 10 10 L & , dust mass of 5 x 10' M Q and 
U-band effective optical depth in stellar birth clouds of 1.6. 
This implies that the typical detected galaxy would be about 
100 times more star-forming, more massive (in terms of stellar 
content) and more dusty, and about 10 times more obscured 
in the optical than the median of the whole UDF sample. 

5.2. Line detections 

We now discuss the possibility of detecting CO and [CII] 
lines in the UDF using frequency scans in all ALMA bands. 
We compute the expected line sensitivities for a t otal on- 
source time of 500 hours on the UDF as in Section 15.11 but 
taking into account the time needed to scan in band in fre- 
quency space. The frequency interval covered in each fre- 
quency setting is Av = 8 GHz. Therefore, the total number 
of settings required to scan a given ALMA band is the to- 
tal frequency range of the band divided by Av (we note that 
the exact setup will depend on the sideband-separations in the 
various bands). The total time spent in each frequency set- 
ting is then the integration time per pointing divided by the 
required number of settings in each band. We assume the 
same mosaicking scheme of the UDF and use the same in- 
tegration time per pointing in each band as for the continuum 
(Section [5.1| ). The total number of frequency settings, time 
per setting and resulting line sensitivity erii nc (computed us- 
ing the ASC and assuming a bandwidth of 300 km s _1 in 
order to resolve the lines in velocity space) are listed in Ta- 
ble [2] In Table [3] we show the predicted number of 5<7 line 
detections in each band given the line sensitivities erii no of Ta- 
ble|2l assumi ng bot h the Milky Way and the M 82 CO SLED, 
as in Section 14.7.11 It is clear that with the integration times 
per frequency setting of Table [2] we predict a relatively low 
number of line detections, compared with the predicted num- 
ber of continuum detections. In bands 2 to 6, we predict a 
minimum of about 20 CO line detections in each band, when 
assuming the Milky Way CO SLED; around 100 detections 
are predicted if the CO is more excited. To get more line de- 
tections, one would need to go deeper, which likely implies 
a compromise with the area of the survey. For example, to 
reach a rms of 10 /iJy in band 6 (approximately 8 times deeper 
than reached in the "default" survey set-up considered so far), 
which would yield at least between 200 and 800 detections 
over the whole UDF, one would need over 400 hours of effec- 
tive integration time for each of th e 81 p ointings. 

As mentioned in S e ction 14.7.11 if we use the 
ISolomon & Vanden Bouj d2005h calibration to convert 
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Figure 10. (a) Distribution flfjlhe velocity-integrated flux of the [CII] line, I\cu] > f° r a " me galaxies in our sample, (b) Distribution of the number of detections 
as a function of line flux, S v (computed assuming a 300 km s _1 line width), expected in the highest frequency bands available for ALMA and PdBI/NOEMA. 



Table 2 

Summary of a possible observational set-up to observed the Hubble UDF with ALMA in full operations with a total on-source time of 500 hours. 
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time per 
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frequency 


primary beam 


number of 


frequency 


pointing b 


frequency setting 


C^cont 


"line 


3<r continnuum 




/ GHz 


/ arcsec 


pointings" 


settings 


/h 


/h 


/ mJy 


/ mJy 


detections' 1 


ALMA 1 


38 


140 


3 


2 


167 


83.5 


3.7 x 10" 4 


7.6 x 10" 3 


36 


ALMA 2 


80 


76 


10 


3 


50 


16.7 


9.9 x 10~ 4 


1.7 x 10~ 2 


155 


ALMA 3 


100 


62 


15 


4 


33 


8.25 


1.5 x 10~ 3 


2.7 x 10" 2 


221 


ALMA 4 


144 


43 


30 


5 


17 


3.40 


2.5 x 10~ 3 


4.2 x 10~ 2 


363 


ALMA 6 


230 


26 


81 


8 


6.2 


0.78 


5.1 x 10~ 3 


8.4 x 10~ 2 


601 


ALMA 7 


345 


18 


169 


13 


3.0 


0.23 


1.3 x 10" 2 


2.3 x lO" 1 


522 


ALMA 8 


430 


14 


278 


15 


1.8 


0.12 


8.0 x 10~ 2 


1.3 x 10° 


136 


ALMA 9 


660 


9.3 


630 


14 


0.8 


0.06 


2.4 x lO" 1 


3.1 x 10° 


84 


ALMA 10 


870 


1.1 


1080 


21 


0.5 


0.03 


7.1 x lO" 1 


9.5 x 10° 


40 



Note. — Sensitivities computed using the ALMA Sensitivity Calculator (50 antennas and default weather conditions). 
a Computed by dividing the total UDF area by the area of the primary beam in each band. 
b Computed by dividing 500 hours by number of pointings in each band. 
c 8 GHz bandwidth. 
d 300 km/s bandwidth. 

c Expected number of continuum detections over the whole UDF field for each band. The expected number of line detections is given in Table|3] 



the infrared luminosities of our galaxies into CO line lu- 
minosities, we would obtain intrinsically brighter CO lines 
by about one order of magnitude. This would make the 
number of predicted CO line detections in the observational 
setup discussed here much higher than the numbers listed in 
Table [3] In the bands with most predicted detections, ALMA 
2 to ALMA 6, the minimum number of detections (when 
assuming the Milky Way SLED) would increase from on 
average 22 to 1 16 in each band, and the maximum number of 
detections (when assuming the M 82 SLED) would increase 
from on average 50 to 242 in each band. 

We use the same method to estimate the number of expected 
[CII] line detections in the highest frequency ALMA bands, 
based on our estimates of [CII] line fluxes from Section l4.7.2l 
We predict a total of 15 [CII] 5cr-detections using the integra- 
tion times per pointing and per frequency setting in each band 
listed in Table|2](Table|4]l. In bands 8, 9 and 10, these numbers 
represent a low fraction (< 0.1 per cent) of the total number 
of galaxies in the observable redshift range because the large 
number of required pointings to cover the whole UDF implies 



a very short integration time per pointing and hence a rela- 
tively high rms. The minimum dust luminosities and star for- 
mation rates of the galaxies that can be detected in [CII] using 
this observational setup are 5 x 10 11 Lq and 10 Mq yr _1 , 
at the very high end of the distribution for the whole sample 
(e.g. FigflTTi. With band 7, the number of detections is higher 
and because the rms is smaller thanks to a higher integration 
time, allowing us to go deeper and detect [CII] emission from 
galaxies with dust luminosities as low as 1.4 x 10 11 Lq. 

6. DISCUSSION & CONCLUSIONS 

In this paper, we have presented empirical predictions of 
(sub-)millimeter continuum and CO/[CII] line fluxes for a 
sample of 13,099 galaxies with redshifts up to z = 5, selected 
in the deepest optical/near-infrared catalog of the Hubble Ul- 
tra Deep Field over 12 arcmin 2 . We have performed a self- 
consistent modelling of the observed optical/near-infrared 
spectral energy distributions of the galaxies, which relies on 
an energy -balance technique, and have allowed us to obtain 
Bayesian estimates of the total infrared luminosity and (sub- 
)millimeter continuum flux densities in several ALMA, JVLA 
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Figure 11. Comparison between the distribution of redshift (a), V-band magnitude (b), stellar mass (c), star formation rate (d), dust luminosity (e) and V-band 
optical depth in stellar birth clouds (f) for our whole sample of galaxies (gray filled histograms), and for the 601 galaxies expected to be detected at the 3tr-level 
in the continuum in ALMA band 6 at 230 GHz, following the observational setup described in Table|2](orange histograms, with numbers given on the right-hand 
y-axis). It is clear from these plots that the detected galaxies will be the most massive, star-forming and dust-rich galaxies in the UDF. 



and PdBI/NOEMA bands. We then combined the constraints 
on total infrared luminosity of our galaxies with empirical 
correlations to obtain estimates of their CO and [CII] line 
luminosities. One advantage of our method is that it allows 
us to derive reliable confidence ranges for the physical pa- 
rameters of the galaxies derived from SED fitting. Using 
fits to the complete ultraviolet to far-infrared SEDs of a sub- 
sample of galaxies in the UDF, we show that, even though 
our confidence ranges for the infrared luminosities and (sub- 
)millimeter continuum fluxes can be as large as one order of 
magnitude, our estimates are reliable, and the large confidence 
ranges reflect the inherent uncertainty in deriving infrared and 
(sub-)millimeter properties of galaxies from optical observa- 
tions. 

Our predictions rely mainly on two assumptions. First, that 
we can reliably predict the total infrared luminosity from the 
attenuated stellar emission, i.e. that our dust attenuation pre- 
scription is correct. We expect this to be true for moderately 
star-forming, optically-selected galaxie s with moderate dust 
optical depths a nd infrared luminos ities ([C hariot & Fal J2000t 
da Cunha et ail I2008t iDaddi et alj 120051 120071: iReddv et alj 
2006). Another argument in favor of the relatively low typical 
optical depths of our galaxies is that, if we assumed infrared- 
to-optical ratios typical of very dust-obscured ULIRGs, we 
would strongly overestimate the extragalactic background 
light at 870/im (Appendix [A). The second assumption is that 
the empirical relations between infrared luminosities and line 
emission, calibrated for low-luminosity galaxies only in the 
local Universe for obvious observational limitations, still hold 
for low-luminosity galaxies out to z ~ 5. With this in mind, 
we have discussed how our results depend on the adopted 
Lir — L' co correlation. To compute the luminosities of CO 
lines with upper level J > 2 from the CO(1-0) line luminos- 
ity, we rely on empirical CO spectral line energy distributions, 
which are a tracer of the excitation state of the gas. Since this 



is a particularly unknown for high-redshift low-luminosity 
galaxies, we discuss predicted CO line luminosities assum- 
ing two extreme CO SLEDs - Milky Way and M82 - which 
should bracket all possible CO excitations in these galaxies 
(assuming star formation is the only excitation source, i.e. 
no AGN). Considering these uncertainties, our predictions are 
the best possible based on the current empirical knowledge of 
star-forming galaxies. 

A possible caveat of the optical selection of our sample 
is that we may be missing very optically-fai nt dust-obscured 
sources, which would not be included in the lCoe et al.l (12006b 
optical catalog but could be bright in the (sub-)millimeter. 
This is most relevant for galaxies with z > 2, where our op- 
tical catalogue becomes incomplete (Fig.Q]). However, based 
on how well our predicted extragalactic background light and 
number counts at 870/im agree with observations, we con- 
clude that such sources are probably not dominant. If very 
dust-obscured sources that are not included in our catalog ex- 
ist in the considered field, our optically-based predictions can 
be considered a lower limit for the number of detections in a 
(sub-)millimeter survey. 

Barring cosmic variance, our predictions should provide a 
statistical representation of an arbitrarily chosen region on the 
sky, and can be used to plan deep field observations in the 
(sub-)millimeter regime. We have illustrated the usefulness of 
these predictions by estimating the expected number of con- 
tinuum and line detections in a 500 hour deep survey of the 
Hubble UDF with ALMA. For example, at 230 GHz (ALMA 
band 6), considering the required number of pointings neces- 
sary to cover the total area of the UDF (given the primary 
beam size), we expect a continuum sensitivity of 5.1 /jjy 
(computed using the ALMA sensitivity calculator, assuming 
a bandwidth of 8 GHz). According to our predictions, this 
would yield a total of ~ 600 3er-detections of galaxies dis- 
tributed uniformly in redshift, with star formation rates down 
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Table 3 

Predicted number of 5 cr CO line detections in the Hubble UDF with ALMA, 
assuming a total on-source integration time of 500 hours, for two differ ent 
cases of CO line SLEDs: Milky Way and M 82 center (SectiongTT). 



ALMA 1 [31.3 - 45 GHz] 

5a-iine = 3.82 X 1CT 2 mJy 



Line 


Milky Way 


M 82 center 


CO(l-O) 


17 


17 


CO(2-l) 


5 


16 




ALMA 2 [67 - 90 GHz] 






5<x linc = 8.55 x 10~ 2 mJy 




Line 


Milky Way 


M 82 center 


CO(l-O) 


21 


21 


CO(2-l) 


7 


16 


CO(3-2) 


1 


38 


CO(4-3) 


2 


14 




ALMA 3 [84 - 116 GHz] 






5(7 Um = 1.34 X 10 _1 mJy 




Line 


Milky Way 


M 82 center 


CO(l-O) 


5 


5 


CO(2-l) 


15 


37 


CO(3-2) 





17 


CO(4-3) 





37 


CO(5^4) 





1 1 




ALMA 4 [125 - 163 GHz] 






5cr linc = 2.10 X lO -1 mJy 




Line 


Milky Way 


M 82 center 


CO(2-l) 


14 


26 


CO(3-2) 


4 


27 


CO(4-3) 





12 


CO(5-4) 





32 


CO(6-5) 





27 


CO(7-6) 





8 




ALMA 6 [211 - 275 GHz] 






5cr linc = 4.22 X lO -1 mJy 




Line 


Milky Way 


M 82 center 


CO(2-l) 








CO(3-2) 


10 


30 


CO(4-3) 


6 


39 


CO(5^) 


3 


43 


CO(6-5) 





18 


CO(7-6) 





7 




ALMA 7 [275 -373 GHz] 






5cr linc = 1.15 X 10° mJy 




Line 


Milky Way 


M 82 center 


CO(3-2) 








CO(4-3) 





22 


CO(5^) 





30 


CO(6-5) 





22 


CO(7-6) 





10 



Table 4 

Predicted number of 5<r [CII] line detections in the Hubble UDF with 
ALMA, assuming a total on-source integration time of 500 hours. 



Band 


Frequency range 


5o"li ne 


number of 5<r 




/ GHz 


/ mJy 


[CII] detections 


ALMA 7 


275 - 373 


1.15 


10 


ALMA 8 


385 - 500 


6.65 





ALMA 9 


602-710 


15.6 


3 


ALMA 10 


787 - 950 


47.3 


2 



to 1 M Q yr _1 and infrared luminosities down to 10 10 L®, We 
also discuss the predicted number of CO line detections when 
performing a scan over the whole frequency range covered 
by each band in steps of 8 GHz (the interval covered in each 
frequency setting). In 500 hours, the line sensitivity reached 
by such a frequency scan of each UDF pointing in band 6 is 
84 /iJy (assuming a 300 km s _1 bandwidth), which would 
yield between 19 and 137 detections, depending on the CO 
excitation. 

The predictions presented in this paper will help to plan 
future (sub-)millimeter line and continuum deep field cam- 
paigns with new cutting-edge interferometers, and will also 
serve as a benchmark of our current empirical knowledge 
of dust continuum emission and CO/[CII] line emission, to 
which these future deep observations can be compared. 
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Figure 12. Ultra violet to millimetre spectral energy distributions of Arp220 (ULIRG), M82 (low-luminosity starburst), and M 100 (spiral galaxy) from 
Silva et al. ( 1998), normalized in the V-band to emphasize their different optical-to-infrared ratios. 
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APPENDIX 

A. COMPARISON WITH OTHER SED TEMPLATES 

For reference, we compare our predicted flux densities with those obtained using three typical local galaxy SEDs: an ultra- 
luminous infrared g alaxy (ULIRG), A rp 220, a low-luminosity starburst, M82, and a typical spiral galaxy, M 100. We use 
template SEDs from Silva et al. ( 1998) which cover the full spectrum from ultraviolet to radio wavelengths. In Fig. [12] we plot 
these SEDs normalized in the F-band to illustrate the large range in optical-to-infrared ratio spanned by these templates. 

We k-correct and re-scale each template to fit the observed optical magnitudes of each galaxy in our sample, and from that we 
predict what the continuum fluxes in the ALMA bands would be if the observed galaxy SED was the same as the template. 

In Fig. [13] we compare the distributi on of the infrared luminosity of the galaxies and the continuum flux in three ALMA bands 
when using our Ida Cunha et all (|2008) SED fits and when using these three templates. These distributions show the range of 
infrared luminosities and (sub-)mm fluxes of our observed galaxies if they all had optical-to-infrared ratios and infrared SEDs 
similar to Arp220, M 82 andM 100, versus our predictions, which make an educated guess by comparing the observed UV/optical 
SED of each galaxy with a library of models with a wide range of optical-to-infrared ratios and infrared SED shapes. It is clear 
that, if we assumed that all observed galaxies had Arp 220-like SEDs, the typical infrared luminosities and predicted continuum 
fluxes would be typically 2 orders of magnitude higher. An M 82 template, with a similar infrared S ED shape but higher op tic al- 
to-infrared ratio, still produces infrared luminosities about one order of magnitude higher than our Ida Cunha et al. (2008J) SED 
fits, but the typical (sub-)mm continuum fluxes are similar to ours, except for part of the sample presenting an excess of fluxes 
compared to our estimates. The typical spiral template, M 100, has the higher optical-to-infrared ratios and cooler typical dust 
temperature, and pre dicts similar, if not slightly lower, typical infrared luminosities and sub-mm continuum fluxes than our fits. 

As in Section 14.51 we can compare the integrated continuum emission over the whole UDF, which, divided by the area of 
the field, allows us to get an estimate of the extragalactic background light (EBL). The predicted continuum fluxes at 345 GHz 
(ALMA 7) from our da Cun h a et al.l ([2008) SED fits yield a EBL of 45.6 Jy deg -2 , which is consistent with observations, 
as discussed in Section |4~51 When using fixed SED templates, we get 389 Jy deg -2 (Arp 220), 52.5 Jy deg~ 2 (M82) and 
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log [ F„ (230 GHz) / Jy ] log [ F, (100 GHz) / Jy ] 

Figure 13. Comparison between our estimates (using the da Cunha et al. 2008 models) and estimates using 3 typical galaxy templates of (a) infrared luminosity; 
(b), (c) and (d) distribution of continuum flux density at 345, 230 and 100 GHz (ALMA bands 7, 6, 3; PdBI/NOEMA bands 4, 3, 1), respectively. 
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redshift 



Figure 14. Observed frequency of the CO lines CO(1-0) to CO(7-6) (from darker blue to lighter blue lines) and [ClI] 158/xm (red line) as a function of redshift. 
The gray and green shaded areas show the frequency ranges sampled by the ALMA, JVLA and PdBI/NOEMA bands considered in this paper. 



9.32 Jy deg 2 (M 100), i.e. in the case of Arp 220 we get too much EBL at 345 GHz by almost a factor of 10, while in the case 
of M 100 we tend to underestimate the EBL at 345 GHz by about a factor of 5. 

B. REDSHIFT RANGES WHERE CO AND [CII] ARE OBSERVABLE IN TYPICAL (SUB-)MM BANDS 

For reference, we show in this appendix the redshift ranges with the CO lines (1-0) to (7-6) and the [CII] line are observable 
in the typical (sub-)mm bands from ALMA, JVLA and PdBI/NOEMA (Fig. [Hand TableEJ. 



Table S 

Redshift ranges where the CO and [CII] lines are observable as a function of band frequency range. 



CO(l-O) CO(2-l) CO(3-2) CO(4-3) CO(5-4) CO(6-5) CO(7-6) [C II] 

115.27 GHz 230.54 GHz 345.80 GHz 461.04 GHz 576.27 GHz 691.47 GHz 806.65 GHz 1900.54 GHz 



Band 



Frequency range 
/GHz 



JVLAK 
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26.5 


3.35 


- 5.40 
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- 12.0 
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1.56 


-2.68 
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-6.37 
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-10.0 


9.24- 


-13.7 
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1.31 


-1.88 


3.61- 
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8.22- 
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67- 


90 


0.28 


- 0.72 


1.56- 


-2.44 
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6.68- 
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C. PREDICTED (SUB-)MM FLUXES FOR EACH GALAXY IN OUR SAMPLE 

In this Appendix, we provide the predicted continuum and line fluxes for each galaxy in our sample. In Table [6] we list the 
predicted continuum fluxes (including confidence ranges) in the various (sub)mm bands. In Table [7] we provide the predicted 
fluxes of the CO lines (using both Milky Way and M82 SLEDs) and of the [CII] line. 



Table 6 

Predicted continuum fluxes and confidence ranges for the first 10 galaxies in our sample. The predicted fluxes are the median of the likelihood distribution computed as described in Section l3~4l and the 

confidence ranges correspond to the 16th - 84th percentile of the likelihood distribution. 



ID R.A. Dec Zphot log(F„/Jy) log(F„/Jy) log(F„/Jy) log(F„/Jy) log(F„/Jy) log(F„/Jy) log(F„/Jy) log(F„/Jy) log(F„/Jy) 

38 GHz 80 GHz 100 GHz 144 GHz 230 GHz 345 GHz 430 GHz 660 GHz 870 GHz 
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-27:49:43.20 


1.029 


-8 


fi2 +0.71 


e no + 0.30 


-b.06±8:|| 


-7 


fi2 +0.41 
DZ -0.63 


-6 


c.3+0.44 
ao — 0.60 


-6 


3O+0.44 


-6 


,0+0.43 
■ lo -0.60 


-5. 


fi o+0.37 
"°-0.65 


-5 


4O+0.32 
• 49 -0.68 



Note. — IDs are the same as in the Coe et al. 12006) catalog. This table is published in its entirety in the electronic edition the journal; a portion is shown here for guidance regarding its form and content. 
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Table 7 

Predicted CO and [CII] line fluxes for the first 10 galaxies in our sample. For the CO lines, we include t he flux prediction using a Milky Way CO SLED (MW) and a M82 center SLED (M82), as described in 

Section gTT] 



ID 


2 phot 


log(S„/Jy) 


log(S„/Jy) 


log(S„/Jy) 


log(S„/Jy) 


log(S„/Jy) 


log(S„/Jy) 


log^/Jy) 


Iog^/Jy) 






CO(1-0) 


CO(2-l) 


CO(3-2) 


CO(4-3) 


CO(5-4) 


CO(6-5) 


CO(7-6) 


[CII] 








MW / M82 


MW / M82 


MW / M82 


MW / M82 


MW / M82 


MW / M82 




1 


0.472 


-4.58 


-4.28/-3.99 


-4.19/ -3.66 


-4.14/ -3.45 


-4.27/ -3.31 


-4.50/ -3.24 


-4.63/ -3.30 


-2.36 


2 


2.703 


-6.64 


-6.34/ -6.04 


-6.25/ -5.72 


-6.20/ -5.50 


-6.32/ -5.37 


-6.56/ -5.30 


-6.69/ -5.36 


-4.50 


3 


1.281 


-7.23 


-6.93/ -6.64 


-6.84/ -6.31 


-6.79/ -6.10 


-6.92/ -5.96 


-7.16/ -5.90 


-7.28/ -5.96 


-5.26 


4 


3.791 


-6.66 


-6.35/ -6.06 


-6.26/ -5.73 


-6.21/ -5.52 


-6.34/ -5.38 


-6.58/ -5.32 


-6.70/ -5.38 


-4.49 


5 


0.452 


-5.83 


-5.52/ -5.23 


-5.44/ -4.91 


-5.38/ -4.69 


-5.51/ -4.55 


-5.75/ -4.49 


-5.88/ -4.55 


-3.78 


6 


2.769 


-6.95 


-6.65/ -6.36 


-6.56/ -6.03 


-6.51/ -5.82 


-6.64/ -5.68 


-6.87/ -5.62 


-7.00/ -5.68 


-4.86 


7 


1.653 


-5.73 


-5.42/ -5.13 


-5.33/ -4.80 


-5.28/ -4.59 


-5.41/ -4.45 


-5.65/ -4.39 


-5.77/ -4.45 


-3.52 


8 


0.636 


-3.94 


-3.64/ -3.35 


-3.55/ -3.02 


-3.50/ -2.81 


-3.63/ -2.67 


-3.86/ -2.61 


-3.99/ -2.66 


-1.60 


9 


0.695 


-6.41 


-6.il/-5.82 


-6.02/ -5.49 


-5.97/ -5.28 


-6.10/ -5.14 


-6.33/ -5.08 


-6.46/ -5.14 


-4.40 


10 


1.029 


-6.27 


-5.96/ -5.67 


-5.88/ -5.35 


-5.82/ -5.52 


-5.95/ -4.99 


-6.19/ -4.93 


-6.31/ -3.98 


-4.18 



Note. — IDs are the same as in the Coe et al. ( 2006) catalog. This table is published in its entirety in the electronic edition the journal; a portion is shown here for guidance regarding its form and content. 
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